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Abstract 

The B — > tttt hadronic matrix element of the chromomagnetic dipole operator Og g 
(gluonic penguin) is calculated using the QCD light-cone sum rule approach. The 
resulting sum rule for (-7riT\08g\B) contains, in addition to the 0(a s ) part induced by 
hard gluon exchanges, a contribution due to soft gluons. We find that in the limit 
nib oo the soft-gluon contribution is suppressed as a second power of l/m^ with 
respect to the leading-order factorizable B — ► tttt amplitude, whereas the hard-gluon 
contribution has only an a s suppression. Nevertheless, at finite mj, soft and hard 
effects of the gluonic penguin in B — > tttt are of the same order. Our result indicates 
that soft contributions are indispensable for an accurate counting of nonfactorizable 
effects in charmless B decays. On the phenomenological side we predict that the 
impact of gluonic penguins on B^ — ► 7r + 7r~ is very small, but is noticeable for B^ — > 
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1 Introduction 



Intensive experimental studies of B — > tttt and other charmless two-body hadronic B decays 
are currently under way at B factories pQ. The major goal is to complement the already 
observed CP- violation in the B — > J / ipK channel by measuring other CP- violating effects 
and quantifying them in terms of the angles a and 7 of the CKM unitarity triangle (for 
a recent comprehensive review, see, e.g. Ref. [2 ). The extraction of these fundamental 
parameters from experimental data on B — ► tttt, Kit, KK demands reliable theoretical 
estimates of hadronic matrix elements of the effective weak Hamiltonian between the initial 
B and the final two- light-meson states. These matrix elements are not yet accessible in 
lattice QCD, and thus other QCD methods of their calculation are actively being developed. 

The factorization ansatz in hadronic two-body B decays has recently been put on a more 
solid ground. It has been shown to 0(a s ), that in the limit m& — > 00 the amplitude for 
B — > tttt factorizes into a product of the B — > tt form factor and the pion decay constant. 
The nonfactorizable effects due to hard gluon exchanges are directly calculated combining 
perturbative QCD with certain nonperturbative inputs. The latter include the B — > tt form 
factor and the light-cone distribution amplitudes (DA) of pions and B meson. According 
to QCD factorization, various soft nonfactorizable contributions to B — > tttt, such as 
exchanges of low virtuality gluons between the "emission" pion and the remaining B — > tt 
part of the process, vanish in the mj, — > 00 limit. An open phenomenological problem for 
QCD factorization is to obtain quantitative estimates of soft nonfactorizable effects at a 
finite, physical b quark mass. Among these effects, the contributions of current-current 
operators in the penguin topology (e.g., the "charming penguins") could be important, as 
was argued in Ref. 

Recently, a new method to calculate the B — > tttt hadronic matrix elements from QCD 
light-cone sum rules (LCSR) HH H] has been suggested by one of us jHj. The main 
advantage of LCSR is that the hadronic matrix elements of heavy-to-light transitions can 
be calculated including simultaneously soft and hard effects. Thus, the B — > tt form factor 
including both soft (end-point) contributions and hard-gluon exchanges was obtained from 
LCSR HOI UH U2j providing the main input for the factorizable B — > tttt amplitude. 
Furthermore, in Ref. [Hj it was shown that with LCSR one achieves a quantitative control 
over l/rrift effects of soft nonfactorizable gluons in B — ► tttt. The soft-gluon effect in 
the matrix element of the current-current operator in the emission topology turns out to 
be of the same size as the hard-gluon nonfactorizable contributions obtained from QCD 
factorization. It is important therefore to investigate other nonfactorizable effects, related 
to the penguin and dipole operators and/or to non-emission topologies for the current- 
current operators. Especially interesting is to compare the size of soft- and hard-gluon 
contributions. 

To continue a systematic study in this direction we apply in this paper the LCSR 
method to the B — > tttt hadronic matrix element of the chromomagnetic dipole operator 
0$ g (gluonic penguin). Apart from being phenomenologically interesting by itself, the 
(7T7r\Osg\B) matrix element provides a very useful study case for the LCSR approach. 
Indeed, as we shall see, the hard-gluon contribution is given by one-loop diagrams and is 
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therefore relatively simple. The soft-gluon effect corresponding to a tree-level diagram is 
also easily calculable. One is therefore able to estimate both hard and soft gluonic-penguin 
contributions using one and the same method and input. 

The paper is organized as follows. In section 2 we introduce the relevant three-point 
correlation function and calculate it using operator-product expansion (OPE) near the 
light-cone and taking into account both hard- (0(a s )) and soft-gluon contributions. In 
section 3, following the procedure described in Ref. |3, we match the result of this calcula- 
tion to the dispersion relations in 7r and B channels combined with quark-hadron duality, 
and obtain the sum rule for the {7nr\Og g \B) matrix element. In section 4 we investigate 
the rrib — > oo limit of the sum rule. Section 5 contains the numerical analysis. Further- 
more, in section 6 our prediction is discussed from the phenomenological point of view 
and compared, in section 7, with the corresponding result of QCD factorization. We con- 
clude in section 8. The appendix contains some useful formulae: the decompositions of the 
vacuum-to-pion matrix elements in terms of the light-cone DA of various twist. 



2 Correlation function 

Following Ref. [H] we start with defining a vacuum-pion correlation function: 

F a ( P ,q,k) = -J *xe**-*>J Ae iH)!, (0 | T{j<£(y)O 8 ,(0)j< fl \x)} \ *-(q)) , (1) 

where the quark currents = u^ a ^d and jg = rribbijsd interpolate n and B mesons, 
respectively, and the quark-gluon chromomagnetic dipole operator 8g relevant for B — > nn 
has the following standard expression: 

O 89 = 0^(1 + 75 )^G;6. (2) 
This operator enters the AB = 1 effective weak Hamiltonian 

+ ... + X t c 8 g(^)0 Sg (n)\ , (3) 

where \ q = V q bV* d , q — u,t are the combinations of CKM parameters, and we have singled 
out the relevant current- current operators: 

1 = (dl>)(«r^) , d 1 = (dT^u)(ur^b) , (4) 

denoting by ellipses all remaining operators which we do not consider in this paper. In 
Eqs. ©-©, Tr(A a A 6 ) = 25 ab , c\ y 2$ g are the Wilson coefficients, /i ~ nib is the normalization 
scale and = 7^(1 — 75). Furthermore, the operator 2 = (uT^u) (dT^b) has been Fierz 
transformed: 2 = Oi/3 + 20\. 

The correlator ([[J is a function of three independent momenta, which are chosen to be 
q, p — k, k, and, for simplicity, p 2 = k 2 = 0. Similar to Ref. jHj, we introduce the unphysical 
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external momentum k flowing into the Og g vertex. At k ^ 0, the momentum p — q in the 
B channel is independent of the total momentum P = p — g — k of the light-quark state 
formed after the b quark decay. That will allow us later to avoid contributions of the 
"parasitic" light-quark states in the dispersion relation in the variable (p — q) 2 - In the final 
sum rule k vanishes in the ground-state B — ► ttti contribution. In Eq.Q the pion is on 
shell, q 2 = m 2 . We will work in the chiral limit and set m n = everywhere, except in the 
enhanced combination fi n = m 2 /(rn u + rrid)- The Lorentz-decomposition of the correlation 
function contains four invariant amplitudes: 

F a = {p- k) a F + q a F\ + k a F 2 + e a( , Xp q f3 p x k p F 3 , (5) 

depending on three kinematical invariants (p — k) 2 , {p — q) 2 and P 2 = (p — q — k) 2 . In 
what follows only the amplitude F is needed. 

In the region of large virtualities: (p — k) 2 , (p — q) 2 , P 2 < 0, \(p — k) 2 \, \ {p — q) 2 \, \P 2 \ ^> 
Aq CD , the correlation function (0) is calculated using light-cone OPE, in a form of the 
sum of short- distance coefficient functions (hard amplitudes) convoluted with the pion DA 
of growing twist and multiplicity. Only a few first terms of this expansion are usually 
retained, the components with higher twist / multiplicity are neglected being suppressed by 
inverse powers of large virtualities. 

The gluon emitted from the vertex Os g in the correlation function contributes in two 
different ways. First, a highly virtual (hard) gluon is absorbed by one of the quark lines at a 
short light-cone separation from the emission point. The corresponding one-loop diagrams 
(see Fig. 1) are then calculated perturbatively, as an 0(a s ) part of the hard amplitude. The 
second possibility is that the emitted gluon together with the quark- ant iquark pair forms 
the quark-antiquark-gluon DA of the pion (see Fig. 2a). We will call these low-virtuality 
gluons soft. Thus, in the correlation function the hard- and soft-gluon contributions belong 
to different terms of the light-cone OPE and can be calculated separately. In what follows, 
we take into account all 0(a s ) hard-gluon effects of twist 2 and 3, and, in addition the effects 
of quark-antiquark-gluon DA of twist 3 and 4 in zeroth order in a s . This approximation 
corresponds to the accuracy of the light-cone expansion adopted for the LCSR calculation 
of the B — > it form factor. 

The higher twist terms corresponding to the soft-gluon contributions are subleading in 
the light-cone OPE, being suppressed by inverse powers of large virtualities. However, in 
the correlation function (Q) the lowest twist contribution is of 0(a s ), therefore both hard 
and soft contributions can be equally important. We may again refer to the LCSR for the 
B — > it form factor where the 0(a s ) twist 2 and the zeroth-order in a s twist 4 contributions 
are at the same level numerically. 

In the following three subsections the calculation of the various contributions to the 
correlation function (£[]) is presented. 

2.1 Hard gluon contribution 

We begin with calculating the leading 0(a s ) contribution to the correlator ((TJ) generated 
by hard gluons. To obtain the relevant diagrams one has to contract the 6-quark fields 
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Figure 1: Diagrams corresponding to the perturbative 0(a s ) contributions to the correlation 
function (Op. Solid, dashed, wavy lines and ovals represent quarks, gluons, external currents 
and pseudoscalar meson DA, respectively. Thick points indicate the gluonic penguin vertex. 



and, in addition, the d- and d— quark fields from the operators and Og g , respectively, 
inserting the free-quark propagators for both contractions. The remaining on-shell u and 
d fields form the quark-antiquark DA of the pion. The decomposition of the relevant 
vacuum-pion matrix elements in terms of DA of the lowest twists 2 and 3 is presented 
in the appendix. Twist 4 components are neglected, because the 0(a s ) twist-4 effects 
are beyond the accuracy of our calculation. The gluon emitted by 8g is absorbed by 
one of the four quark lines leading to the four diagrams shown in Fig. 1. Note that we 
disregard the possibility to contract the d quark from Os g with the d quark from j\ . That 
type of contraction leads to "penguin-annihilation" diagrams where the heavy-light loop is 
connected by a single gluon with the pion part of the correlation function. An additional 
gluon has to be added in this case to obey color conservation, making this contribution 
either 0(a%), with both gluons hard, or 0(a s ), with one hard gluon and one soft gluon, 
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the latter entering the pion DA. Both effects are beyond the adopted approximation for 
the correlation function and presumably very small. 

The one-loop diagrams are calculated employing dimensional regularization. Among 
other tools we used the FORM program ^H]- The invariant amplitude F is obtained 
by taking the coefficient at (p — k) a . The result for F is ultraviolet divergent; however, 
this divergence does not play any role in the resulting sum rule because the 1/e terms 
vanish after Borel transformations. On the other hand, we find that the amplitude F does 
not contain infrared-collinear divergences. This is in accordance with the fact that the 
perturbative expansion for the correlation function starts, in twists 2 and 3, with 0(a s ). 
If one takes into account an additional perturbative gluon correction in the correlation 
function, the resulting two-loop diagrams will yield infrared-collinear divergences which 
should be absorbed by the evolution of the pion DA, as in the case of the 0(a s ) correction 
to the LCSR for the B -> vr form factor [Till HQ. 

According to the procedure of the sum rule derivation explained in Ref. [S] we have 
to retain in the final answer for the amplitude F only those terms which, in the limit of 
large \P 2 \ ~ m 2 B , have a nonvanishing double imaginary part in the variables (p — k) 2 and 
(p — q) 2 , taken in the duality region < (p — k) 2 < sj, m 2 < (p — q) 2 < s$ . We find 
that only two diagrams, in Fig. la,b, contain nonvanishing terms. Retaining only these 
terms we obtain the following analytical expression for the hard-gluon contribution to the 
amplitude F: 

F hard ( Sl , s 2 , P 2 ) = F tw2 ( Sl , s 2 , P 2 ) + F tw3 ( Sl , s 2 , P 2 ) , (6) 
where the twist 2 and 3 parts are, respectively 
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(8) 



In the above, we denote s± = (p — k) 2 , s 2 = (p — q) 2 and u = 1 — u for brevity. The 
terms proportional to ln(— S\/ p 2 ) and ln((m^ — s 2 )/ml) originate from the diagrams in 
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Fig. la and lb, respectively. Note that the logarithmic dependence on the dimensional 
regularization scale \x is unimportant because it vanishes after Borel transformation in the 
variable s\. 



2.2 Soft gluon contribution 

In order to reach the same accuracy in the light-cone expansion of the correlation function 
(JIJ as in the calculation of the B —>■ n form factor one has to take into account the 
contributions of quark- ant iquark-gluon three-particle DA. Being a regular part of the light- 
cone expansion, these contributions correspond to a nonperturbative effect of zeroth order 
in a s which is qualitatively different from hard-gluon exchanges. The on-shell gluon field 
emitted by the penguin operator is directly absorbed in the three-particle pion DA yielding 
a single tree-level diagram shown in Fig. 2a. The corresponding expression is obtained by 
contracting b and d quark fields and using the decomposition of the (0\u(y)G IJiU (0)d(x)\7r) 
matrix element in terms of twist 3 and 4 quark-antiquark-gluon DA given in the appendix. 
We find that the contribution of twist 3 DA vanishes in the amplitude F and the answer 
contains only twist 4 DA : 



F soft, ^ p2x = <U f VOLi 

1 h 2 ' ' 8vr2 J (ml - (1 - ai )s 2 )(-a 2 pi - (1 - a 2 ) Sl ) 
x{2{ Sl -s 2 - P 2 )[ V± { ai ) + <p±(ai))] + (si + s 2 - P 2 )b||(ai) + £||(oi)]} , (9) 



where a, = a 2 , a^} and T>oti = da\da 2 da^ (1 — ot\ — a 2 — a 3 ) 



2.3 Factorizable 4-quark contribution 

Continuing the light-cone expansion of the correlation function (£[]) one encounters also 
four-particle Fock states. In particular, four-quark (qqqq) states emerge when b quark and 
hard gluon propagate at short light-cone separations and four light-quark fields are on- 
shell, yielding 0(a s ) contribution. One of the corresponding diagrams is shown in Fig. 2b. 
The Fock states with > 3 multiplicity contribute in two different ways. First, various 
multiparticle DA of higher twist can be formed. Usually, in LCSR applications their con- 
tributions are expected to be strongly suppressed (see, e.g. a more detailed discussion in 
Ref. [H] ) and are therefore neglected, also in our analysis. Second, the configurations of 
multiple quark and gluon fields sandwiched between the vacuum and pion state, can in 
general be splitted into a product of two operators of lower multiplicity, one of them being 
vacuum averaged, and the other one forming a vacuum-pion matrix element. In order to 
avoid double counting, these factorizable contributions have to be subtracted from genuine 
multiparticle DA. In our case, the four-quark configurations emerging from the correlator 
(0) contain a factorizable piece, a product of vacuum-averaged quark-antiquark fields and a 
low-twist (2,3) quark-antiquark DA. We will take the four-quark factorizable contributions 
into account approximating the vacuum average of quark and antiquark fields with the 
standard quark condensate of dimension 3. Combined with the twist 2 and 3 DA, these 
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Figure 2: Diagrams corresponding to the soft-gluon (a) and factorizable J^-quark (b) con- 
tributions to the correlation function (QJ). 

contributions are effectively of twist 5 and 6 respectively, and are analogous to the factoriz- 
able twist 6 terms taken into account in the LCSR for the pion electromagnetic form factor 
derived in Ref . ^3] . Other factorizable multiparticle contributions, such as qqGG or qqqqG 
which can be parameterized in terms of dimension 4 (gluon) and 5 (quark-gluon) conden- 
sates are neglected. The purpose of taking into account the quark-condensate contribution 
is twofold. First, it is enhanced by the numerically large parameter fi n , therefore the size 
of the quark-condensate contribution to LCSR serves as an upper limit for all neglected 
multiparticle contributions. Second, the QCD factorization answer for the gluonic penguin 
amplitude contains a certain counterpart of the quark-antiquark condensate term, with 
which our result will be compared. 

Turning to the actual calculation, we find that only one diagram (Fig. 2b) provides 
nonvanishing contributions to the sum rule. The short- distance part of this diagram con- 
sists of b- and u- quark propagators. The momenta of the vacuum averaged <i<i-quarks 
are neglected in the adopted condensate approximation and the vacuum averaging is done 
in a usual way: (0\d^dij\0) = S 1 '- 1 '5^{qq) /12, where (qq) is the quark condensate density, 
q = u,d, and i,j, and £,u are the colour and spinor indices, respectively. The diagram in 
Fig. 2b yields: 



IziT J Si{ — ur — usijyrrifo — US2) 
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3u(s\ — s 2 — P 2 ) — 2si (p p (u) 

{s 1 - P 2 )(4si - 3us 2 ) - 5u{st - P 2 ) 2 + 6sis 2 
—uP 2 — USi 




(10) 



ml - us 2 



where only the terms which contribute to the sum rule are retained. Adding together the 
contributions presented in Eqs. (J%])- pO|) . we obtain our final result for the amplitude F: 



which will be used to derive LCSR. 

3 Dispersion relation and sum rule 

The procedure of deriving the sum rule for the B — > tttt matrix element from the correlation 
function (JIJ is independent of the operator Og g . Therefore we refer to the paper jH] where 
this procedure was explained in detail for a generic case. For completeness, let us summarize 
the main steps of this derivation: 

1. The calculated correlation function (the amplitude F in Eq. ([11))) is matched to the 
dispersion relation in the pion channel, in the variable s\ = (p — k) 2 , at large s\ < 0. 

2. The Borel transformation, si — > M\, is performed and quark-hadron duality is used 
to approximate the contribution of excited pseudoscalar and axial states (with an 
effective threshold sj) yielding, at large spacelike P 2 ~ — ml, an expression for the 
matrix element of the operator product Os g ]B between pion states. 

3. The matrix element (7r|0 8s js|7r) is expanded in powers of the small ratio Sq/(— P 2 ) 
retaining the leading term of this expansion. In the resulting expression the analytic 
continuation from large spacelike P 2 to large timelike P 2 = m 2 B is performed. Impor- 
tantly, in the adopted approximation for the correlation function the matrix element 
(7r7r|O8t,js|0) \ P 2 =m 2 has no complex phase, as directly follows from the calculation 
of the imaginary part in S\ of the relevant diagrams in Figs. 1,2. At < Si < Sq 
the resulting expression remains real at P 2 — > m 2 B . An imaginary part will appear if 
one adds hard gluon exchanges to the diagrams in Figs. 1,2, proceeding beyond the 
adopted order in a s . Physically, the imaginary part of (mi\OsgjB\0)\p 2 =m 2 B should be 
identified, within the accuracy of the quark-hadron duality approximation, with the 
phase of the final-state rescattering of two pions. Our result indicates the smallness 
of this phase. 

4. The expression obtained for (vr7r|O 8ff jB|0) is then equated to the dispersion relation 
in the S-channel, in the variable s 2 = (p — q) 2 - The auxiliary momentum k vanishes 



F(s u s 2 ,P 2 ) 



F hard { Sl , s 2 , P 2 ) + F soft { Sl , s 2 , P 2 ) + F^\ Sl , s 2 , P 2 ) 



(11) 
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in the ground state S-meson contribution. The latter contains the hadronic matrix 
element of Og g between B and irir states multiplied by the 5-meson decay constant 

5. The second Borel transformation, S2 — > Mf , is applied to the resulting relation and 
quark-hadron duality is employed again, with the corresponding threshold s B , to 
approximate the contribution of the excited B states. Finally, the desired sum rule 
for the {7T7r\Os g \B) matrix element is obtained. 

For convenience, we present the resulting LCSR in a form of the sum of three separate 
contributions: 

A<P»\B° d - 7T+7T-) = (n-(p)n+(-q)\O Sg \B d (p - q)) = A^J + A^tf + A^f , (12) 
where 
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(14) 



.(15) 



In the above, f' a {u) = d<p a (u)/du, u B = m\/s B and f B is the B-meson decay constant 
defined as m^^qi^^B) = m 2 B f B . Note that the first bracket in Eqs. ()13j) and (jUj) is 
approximately equal to f n if one uses the SVZ sum rule [To] : 



(0\a s /nG%G a ^\0) 176 
12M 2 + 81M 4 



ira s (qq)'' 



(16) 



retaining only the leading-order quark-loop term. The quark-condensate contribution to 
Eq. (|16j) of the form (m u + md)(qq)/M 2 is absent in the chiral limit, whereas the per- 
turbative correction and the gluon and 4-quark condensate contributions are beyond our 
approximation. 
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The second and third lines in Eq. (|13|) contain the contributions of the diagrams in 
Fig. la and lb, respectively. In the latter contribution only twist 3 DA remain, the twist 
2 term proportional to the small ratio s^jm\ is neglected. For the same reason, in the 
soft-gluon contribution (jUj) only two DA, (p± and (p± contribute. 

To minimize theoretical uncertainties we will calculate the ratio of the gluonic penguin 
matrix element (|T2"|) to the factorizable B — > nir hadronic amplitude. The latter, as shown 
in Ref. jHj, simply coincides with the matrix element of the operator 0\ in the emission 
topology: 

4 0l) (^° - vr+TT-) = (<K-{p)<K + {-q)\O x \Bl{p - q)) E = im'M&iO) , (17) 

if one neglects the 0(sQ/m 2 B ) terms and the two-gluon nonfactorizable exchanges which 
are beyond the adopted approximation. In the above, f Bn (0) is the LCSR result for the 
B — > it form factor at zero momentum transfer. For completeness we write down this sum 
rule P HUl [HJ: 



i 



2m%f B J u \ m b 



ip a {u) uip' a (u 
wp P W + — — 



+ O(qqtwA) + O(qqG) + 0(a s 




where, for brevity, the subleading contributions of the twist-4 quark- ant iquark DA, twist- 
3,4 quark-antiquark-gluon DA and the 0(a s ) correction to the twist 2 term are not shown 
explicitly but taken into account in the numerical analysis. The complete expression for 
LCSR (fTHj) can be found, e.g. in the review jTH]- The recently calculated [12] 0(a s ) 
correction to the twist-3 part of the sum rule (fTSj) is very small and we neglect it here. 



4 Heavy quark mass limit 

The sum rule (|12|) is derived in full QCD, with the finite 6-quark mass. To obtain the 
heavy-quark mass limit of the matrix element A^° S9 \B^ — > 7r + 7r~), we substitute in LCSR 
the expansions of all m^-dependent quantities: 

m B = m b + A , Sq = ml + 2m b u , M 2 2 = 2m b T , f B = f B / y/mj,, (19) 

where A, u Q , r, and f B are the parameters independent of the heavy-mass scale. One 
also has to take into account that, at m b — > oo, u° B —>■ 1 — LO§jm b and only the end-point 
regions of the momentum fraction u contribute in the integrals in Eqs. p3jl -(|15 jl . Using a 
convenient integration variable p = um b , and taking into account the end-point behavior 
of DA, the following substitutions have to be done in these integrals, to the leading order 
in 0{l/m b ) : 

1 2u) 

J du= (l/m 6 ) J dp, (firA u ) = -{p/ m b)^A l ) ' = x > 

o 

<p ± (l - u, 0, u) =(p±{l- u, 0, u) ~ -(p/m b )ip' ± (0, 0, 1) , (20) 
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where <p'±(0, 0, 1) = [d<p±(l — u, 0, u)/du\ u =\. With these substitutions it is easy to repro- 
duce, e.g. the leading l/m?J 2 behavior [7j of the B — > tx form factor from the sum rule 
(|I8j) at m b — > oo: 



mj-too 
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+ 0(l/mf) 



(21) 



where [T7j 



/tt „A/t 



2(^0 



2/ 



(ip e 



p^(l) + /U^(l)--^(l) 



6 



f22l 



is an effective, m^-independent form factor. Consequently, the heavy-quark mass limit for 
the factorizable amplitude obtained from Eq. (|17|) is 
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Applying Eqs. ()19j) and ()20|) to the sum rule ()12j) one obtains the limiting behavior for the 
separate contributions: 
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(24) 
(25) 
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20JQ 
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(26) 



where only the leading in 1/mb term in each of the three contributions is retained. We 
see that the hard-gluon contribution (J23|) has only 0(a s ) suppression with respect to the 
factorizable amplitude (|23)1 . Note that in Eq. (J21|) only the diagram in Fig. la contributes, 
because the twist 3 contribution of the Fig. lb diagram in Eq. f!13|) has an additional 
1/mb suppression. Importantly, the soft gluon contribution ([25)1. where the pion DA is 
proportional to the dimensionful parameter 5%, is suppressed by two powers of 1/mb ■ 
Finally, the quark-condensate contribution (^Bj) is both 0(a s ) and 0(l/m&) suppressed. 
Thus, in the heavy-quark mass limit the sum rule (|12j) behaves similar to the LCSR for 
the B —>■ tx form factor: the twist 2 and 3 DA survive in the asymptotic limit, whereas the 
higher- twist contributions are suppressed by powers of l/m&. 
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Furthermore, we use the relation between the twist 3 DA: f f a {u) = 6(1 — 2u)tp p (u), 
which follows [21] from QCD equations of motion in the approximation where the quark- 
ant iquark-gluon DA are neglected. Hence, at u — > 1, ^ p (l) = — ^(l)/6, and the expression 
in the second line in Eq. (J2U) (and the same in Eq. (|2l)])) is simply equal to /^.(O) given by 
Eq. (j22j). As a result, at — > oo, the gluonic penguin amplitude is approximately equal 
to the heavy-quark limit (J2*3*j) of the factorizable amplitude multiplied by an 0(a s ) factor: 



.4 



(Os s ) 
hard 



m-i, — >oo 



mj, — >oo 



2tt 







2tt 



(27) 



mi- >oo 



where the approximate relation (|T6*j) is used. Note that the above factorization will be 
violated by 0(1/to&) corrections to Eq. (j24j). The latter corrections can be easily obtained 
expanding Eq. (|13|) to the next-to-leading order in the inverse heavy quark mass. Interest- 
ingly however, the quark-condensate contribution (J26|) . being 1/rrib suppressed with respect 
to Eq. ()24|). nevertheless reveals the same factorization property: 



.4 
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yMbfnfBnW = — a I A 



7T + 7T 



f28) 



no- 



where the PCAC relation fi n = —2(qq)/f% is used to replace the quark-condensate density 
by fi n . Naturally, the soft-gluon contribution ()25j) also violates the factorization relation 



(o 8g ) 

hard 



it has no a. 



m b - 



(|27j) . Although this contribution is of order 1/ml with respect to A 

suppression and therefore should be considered separately from the 0(a s /ml) corrections 
stemming from the heavy-mass expansion of hard-gluon and quark-condensate terms in 
LCSR beyond 0(l/mf,). A direct evaluation of the latter corrections from Eqs. (|13|) and 
(115)) will, however, be incomplete because, in deriving LCSR, we have already neglected 
small (calculable) sljrr? B terms of the same order. 



5 Numerical results 

Let us return to the sum rule ()12|) at the finite 6-quark mass and use it to obtain a numerical 
estimate for the ratio 

r (o 8s ) ( £0 ^ 7r+7r - ) = A (o 8g )^o _^ 7l + 7l -)/A { ^ 1 \B Q d -> 7T+7T-) , (29) 

which determines (up to the known Wilson coefficient cs g ) the gluonic-penguin correction 
to the factorizable B — > titt decay amplitude. The ratio (j2Hj) has an important advantage 
of being less sensitive to the input parameters than the individual matrix elements A^ ^ 
and A^ 1 ^. Moreover, r^° s ^ is independent of since the latter cancels in the ratio of 
the sum rules f)12|) and (fTHJl . 
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The parameters for the pion channel are f n = 132 MeV and Sq = 0-7 GeV 2 , the 
latter determined from the two-point SVZ sum rule for /„.. The corresponding Borel 
parameter interval is M 2 = 0.5-1.5 GeV 2 , accommodating also the range used in the LCSR 
for the pion form factor j!4j . 

The inputs for the B channel are taken from the earlier analysis of LCSR (fT%|) for the 
B — > tt form factor in Refs. [HIGH]. In particular, to be consistent with the latter sum rule, 
we identify m& with the one- loop pole mass and adopt the interval m& = 4.7 ± 0.1 GeV 
which is in accordance with the most recent average of the MS mass rhb(rhb) = 4.24 ± 0.11 
GeV [IB] obtained from various model-independent determinations (including QCD sum 
rules for the T system). Fixing the m& interval one determines the corresponding duality 
threshold = 35 =f 2 GeV 2 from the QCD sum rule for fs taken with 0(a s ) accuracy. 
Furthermore we take the range Mf = 8-12 GeV 2 for the second Borel parameter and 

fib = \Jm 2 B — m 2 ~ 2.4 GeV for the related scale at which the pion DA are normalized. 
The same scale is adopted for the a s normalization. For the latter we also use the two-loop 
running with A< 4 > = 280 MeV. However, one should have in mind that in the new sum rule 
(112}) two independent Borel parameters M\ and M 2 are present, and the scale M\ ~ 1 GeV 
is independent of m^,. To investigate the sensitivity to the normalization scale, we will vary 
it in a rather wide interval, between /ib/2 and 2 fib- For the quark condensate density we 
take (qq) (lGeV) = (-240 ± 10 MeV) 3 , or, equivalently, /i^(lGeV) = 1.59 ± 0.2 GeV. The 
normalization parameter of the twist 4 DA £ 2 (lGeV) = 0.17 ± 0.05 GeV 2 is determined 
from the two-point QCD sum rules JTHJ (see also Ref. |2"U]). 

For our initial numerical illustration all pion light-cone DA are taken in the asymptotic 
form, that is, the coefficients a2, /W, u 3n , e n are put to zero in the expressions for DA 
presented in the appendix. The numerical results for the ratio r^° 89 \ together with the 
hard-, soft-gluon and quark-condensate contributions, calculated at the central values of 
the input parameters and for the asymptotic form of the pion DA are plotted in Fig. 3 
as a function of Borel variables. We observe a nice stability in both M\,M2- The most 
important feature of our calculation which is clearly seen in Fig. 3 is that the soft-gluon part 
of the gluonic penguin amplitude, while being suppressed by two powers of mj, numerically 
amounts to about 50% of the hard-gluon contribution and has an opposite sign. A relatively 
large magnitude of the soft-gluon term ()14|) can be traced back to the large numerical factor 
of 20 originating from the contributing twist-4 DA and compensating the suppression 
factor <5 2 /m 2 . The quark-condensate contribution, on the other hand, is sub leading also 
numerically, not larger than 30% of the hard-gluon contribution. 

To take into account deviations from the asymptotic form, we have repeated our cal- 
culation with the nonasymptotic part of the pion DA taken to the next-to-leading order 
in the conformal spin [21]. In this case we adopt in the pion twist-2 DA <22(lGeV) = 
0.24 ± 0.14 ± 0.08 obtained [20] by fitting the LCSR for the pion form factor to the experi- 
mental data. The relevant parameters in the twist 3,4 DA have been estimated at the scale 
of ~ 1 GeV using two-point QCD sum rules: f 3w = 0.0035, lu 37T = —2.88 [22] and e n = 0.5 
[2*Tj . These estimates have an uncertainty of ±30%. The corresponding LO anomalous 
dimensions can be found in the appendix. 
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The difference between the values of r^° 89 ^ obtained with nonasymptotic and asymptotic 
DA is inessential for the hard-gluon and quark-condensate contributions. On the contrary, 
in the soft-gluon term, the effect of the nonasymptotic part determined by the parameter 
e n is substantial. The reason is that in the soft-gluon term described by the diagram 
in Fig. 2a the gluon carries the dominant fraction = u (v,q < u < 1) of the pion 
momentum, the <i-quark carries a.\ — 1 — u, whereas the u-quark momentum fraction is 
restricted to small values: ct2 ~ s^/m^. The integration over the variables on in this 
subspace of the integration region YL% a % = 1 is strongly influenced by the presence of the 
second polynomial in the conformal expansion of the twist 4 DA <p±(ai) and (p±_(ot.i). As 
a result, the soft-gluon contribution, being —50% of the hard-gluon part in the case of 
asymptotic DA, drops to (-10%) -(+20%) at e w = 0.35-0.65. 

Finally, we get the following numerical estimate for the ratio of the gluonic-penguin 
and factorizable amplitudes including nonasymptotic effects: 

r (Os g ) (£0 _^ = 0.035 ± 0.015. (30) 

In obtaining the above range, the uncertainties caused by the variation of all parameters 
within their allowed intervals are added linearly. The uncertainty in Eq. (fHUJ) is mainly due 
to the sensitivity to the nonperturbative parameters e n , (qq) and 5%, whereas the choice of 
the 6-quark mass, Borel parameters and has a smaller impact. There is an additional 
uncertainty of approximately ±20% from varying the normalization scale in the limits 
/if>/2 -2/i5. The LCSR prediction for the B —>■ it form factor with the same input and the 
same conservative procedure of estimating the uncertainty is f^ n = 0.28 ± 0.06. 

One can also investigate the quality of the m& — > oo limit for the ratio r^° 89 \ Dividing 
Eq. (|2~3*|) by Eq. ()27j) we find that, within uncertainties of our input, r^° s ^ varies between 
0.6 and 1.1 of its heavy quark limit. Future improvements of nonperturbative parame- 
ters determining the vacuum-pion matrix elements of twist 3,4 will allow to decrease the 
uncertainty of our numerical estimates, in particular, of the subleading in l/m& terms. 



6 Gluonic penguins and the decay amplitude 



We are now in a position to assess the role of the gluonic-penguin matrix element in 
B° — ► 7r + 7r~. The decay amplitude can be represented in the following form: 



A(B 



7T+7T-) EE (7l + 7l-\H W \B° d ) 



7T 7T 



+ ... + V 8 ,( / i)r( o ^( J B d ^7r + 7r-)| 



(31) 



where we denote by ellipses all terms that are not relevant for the present discussion: the 
hadronic matrix elements of 0\ t 2 i n the annihilation, penguin and penguin-annihilation 
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topologies; the contributions generated by current-current operators with c-quarks ( "charm- 
ing penguins") and, finally, all effects of quark- and electroweak-penguin operators (for a 
model-independent classification of all contributions and topologies see Ref. [22] )• The sec- 
ond line in Eq. (f3Tj) contains the factorizable amplitude and the nonfactorizable correction 
due to the operator 0\. The latter is parametrized by the ratio 

rlPW - tt+O = Af\B« - n+n-)/A^\B° - n^) , (32) 

where A^* 1 ^ is the B d ~ * 7r+7r_ matrix element of 0\ in the emission topology. Finally, in 
the third line in Eq. (pTTj) we added the gluonic penguin correction calculated above. 

The matrix element of 0\ contains the hard 0(a s ) and soft 0(l/m&) contributions, so 
that 

r jfO(fl° _> vr+vr-) = r f > hard \B° d -> tt+tt-) + A ^(^° ~> ^ (33) 

"'/; 

where the soft contribution was calculated from LCSR in Ref. |Sj with an estimate Xe{B^ —>■ 

7r + 7r~) = 0.1 ± 0.05 GeV. The hard nonfactorizable correction r ^ 1,haTd ^ i s also accessible 
with LCSR but demands two-loop calculation. Following Ref. [Hj we simply use the QCD 

factorization estimate: r^ 1,hard \B^ — > 7r + 7r~) = a s CFF/(8irN c ), where the expression for 
F is given in Ref. [3]. Since we are only interested in an order-of-magnitude estimate, we 

neglect the twist-3 part of r ^? 1,/iard ) which in QCD factorization [3] has to be regularized 
introducing some additional parameter. For consistency we use the same normalization 
scale fib and twist- 2 pion DA as in LCSR. We obtain: 

r ( E l] (B° d -> 7T+7T-) = -[(0.002-0.016) + 0.045i] + [0.01-0.03] , (34) 

where the first and second brackets contain the hard-gluon (QCD factorization) and soft- 
gluon (LCSR) contributions. The relatively large uncertainty of the hard-gluon part is due 
to the variation of 02(1 GeV) in the adopted range 0.02-0.46. Substituting the estimates 
(|30|) and (|34|) in Eq. ()31|) and using also the LCSR estimate for the B — > tt form factor we 
obtain: 

A(B° d -> 7T+7T-) = £j=U ^1(0.28 ± 0.05)|A u (l.03 

+ [0.001 -0.009 + 0.025i] - [0.005-0.015]) + ... - A t [0.003 -0.008] \ , (35) 



where, for consistency, the Wilson coefficients [21] are taken at the scale fib- We use 
Ci(fib) = 1-124 and c^fib) — —0.272 in NLO and in the NDR scheme and, correspondingly, 
cs g {fib) = —0.166 in LO. In Eq. (j33|) the first (second) bracket in the term proportional 

to X u is due to the hard (soft) parts of r% . Note that some of the uncertainties in the 
B — > 7r form factor and in the nonfactorizable contributions are correlated. 

Both nonfactorizable corrections in Eq. (J35j) turn out to be very small, not exceeding 
a one-percent level. However, before one can claim that B d — > 7i + tt~ has a predominantly 
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factorizable amplitude, it is necessary to calculate all remaining contributions, having in 
mind that the sum of all nonfactorizable contributions may accumulate into a noticeable 
effect. As both analyses of Ref. jH| and this paper show, in this calculation the soft-gluon 
effects are indispensable. 

Note that the situation with the colour-suppressed B% — > channel is quite different. 
Here the factorizable part is proportional to the small coefficient c 2 (//&) + Ci(jU&)/3 ~ 0.1, 
therefore the contribution of the 0\ matrix element multiplied by 2c x is equally important. 
Also the gluonic-penguin correction in this channel can reach a few % of the total amplitude. 
One may expect that in some other penguin-dominated channels of charmless B decays, 
such as B — > Kn, the effects of gluonic penguins may even be more important. 

7 Comparison with QCD factorization 

It is interesting to compare our estimate of the hadronic matrix element (7m\Osg\B) ex- 
pressed via the ratio r^° 8 ^ with the same ratio derived from QCD factorization. We can 
only compare the 0(a s ) contributions in our calculation because, in the QCD factorization 
scheme, the soft-gluon effects have not yet been investigated quantitatively. 

To simplify the comparison, we set to zero all Wilson coefficients in Hy/ except c§ r In 
this case, the hard-scattering kernel in 0(a s ) calculated within QCD factorization approach 
contains a single diagram with a hard gluon emitted by the 0$ g vertex and converted into a 
quark- ant iquark pair. This diagram (the sixth one in Fig. 2 of the second paper in Ref. [3]) 
is simply reproduced from the one in Fig. la of this paper if one puts on shell the quarks 
emitted by the 5-meson and pion currents. With the same procedure of putting on shell the 
quark lines applied to the Fig. lb diagram, one gets a diagram which describes a hard-gluon 
exchange between the gluonic-penguin vertex and the spectator quark in the 5-meson, with 
a subsequent hadronization of the light-quark pair into a pion pair. The long-distance part 
of this process can be described by a two-pion DA, an object different from the one-pion DA 
used before. In the framework of QCD factorization this particular diagram has not been 
taken into account. A convincing reason for that is that the corresponding contribution is 
l/m b suppressed, which is also confirmed by the heavy-mass expansion of LCSR discussed 
above. We conclude, that, apart from soft-gluon effects, at finite there are certain 
hard-gluon effects in LCSR, such as the contribution of the Fig. lb diagram, which do not 
appear in the QCD factorization. 

After this qualitative discussion we turn to the quantitative comparison with the gluonic 
penguin contribution in the QCD factorization. In our toy scenario with only cg 9 7^ it 
is easy to extract this contribution from the expression for the effective coefficients 
obtained in Refs. [3] for B — > kit. The only effective coefficients where cg g contributes are 
04 and ci6, so that the ratio r^° 89 ^ in the QCD factorization scheme is : 

r { 83 '{B d -> 7T+7T ) QC Dfact. = ^TW / dU l + ' 36 

2nN c \Jo 1 — u m b J 
together with the 0(n n /mb) correction which is the only 1/m?, effect retained in QCD 
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factorization [3] because of the large numerical value of fx n . At m& — > oo the above ratio, 
for the asymptotic DA ip n = 6m (1 — u), coincides with the heavy-quark mass limit (|27|) 
obtained from LCSR, if in the latter the approximation (|T6|) is used. Thus, the LCSR 
prediction for the gluonic-penguin matrix element reproduces the QCD factorization result 
for the same matrix element in the limit of infinitely heavy b quark. 

The two methods, however, substantially differ in subleading terms. This difference 
manifests itself already at the 0(a s /m&) level. In addition to the penguin-spectator ex- 
change effect corresponding to the diagram in Fig. 2b and discussed above, there is a factor 
of 2 difference between the second, subleading term in Eq. (J3fi|) and the heavy quark mass 
limit (jZEj) of the quark-condensate contribution in LCSR. Furthermore, the fact that the 
soft gluon corrections in LCSR are suppressed by two powers of l/m b is in a general agree- 
ment with the expectations of QCD factorization. However, as we already noted, this effect 
which is absent in the QCD factorization prediction, turns out to be essential at finite mj. 

To complete our comparison with other methods let us parenthetically note that LCSR 
predictions do not support the PQCD approach (2H| to charmless B decays. In particular, 
according to LCSR and contrary to PQCD, the soft contributions dominate in the B — > ix 
form factor. 

8 Conclusion 

In this paper, we have presented the first calculation of the hadronic B — ► titt matrix ele- 
ment of the gluonic penguin operator 0$ g using QCD LCSR. This matrix element is a very 
suitable study object for the sum rule approach because both hard- and soft-gluon contri- 
butions are calculable within one procedure using the same input. The results obtained 
for (-K + 'n-^\Osg\B^) from LCSR clearly indicate that at finite m;, the soft-gluon and other 
1/rrib suppressed effects are important motivating further investigation in this direction, 
in particular, calculating the penguin matrix elements of current-current operators. Im- 
portantly, in the heavy-quark mass limit the LCSR prediction is consistent with the QCD 
factorization. Similar to the nonfactorizable emission, the gluonic-penguin contributions 
to B — > mi are suppressed either by a s or by powers of 1/m.f,. On the phenomenological 
side, our result clearly indicates that gluonic penguins play an insignificant role in the 
B% — > 7r + 7r~ decay amplitude but become noticeable in B® — > 7r°7r°. Our prediction can 
also be used in the future analyses of penguin-dominated B — > Ktt modes. 
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Appendix: Vacuum-pion matrix elements and distri- 
bution amplitudes 

The following convenient decomposition of the vacuum-pion matrix element near the light- 
cone is used in the calculation of the correlation function (fT|): 



<O|«t(a; a )40ri)|7r-(g)) 



i5 ij 



12 



fir / du e 



-iuqxi — iuqX2 



+ (lb)^-K i Pp(u) - - (^75)^^/3 (Xi - X 2 ) T fi 7T iPa(u) 



twist > 4, (Al) 
c, x 2 = 0), i,j = 1,2,3 



where the quark fields are taken near the light-cone (xi = 
and £, u — 1 -4 are the quark color and spinor indices, respectively. Only the twist-2 DA 
<£V and the twist 3 DA tp p and <p a are retained. The familiar definitions of these DA are 
obtained by multiplying both parts of this equation by the corresponding combinations of 7 
matrices and taking Dirac and color traces. For the twist-2 pion DA, in the next-to-leading 
accuracy in the conformal spin j2Ij adopted in this paper, only the first nonasymptotic 
term is retained in the expansion in Gegenbauer polynomials C 2 ^ 2 : 

<p n (u, jj) = Quu 1 + 0,2(1*1) C ^ 2 (u — u 



The LO scale-dependence of the nonasymptotic part is determined by 



(0) in 

a 2 (/i 2 ) = [Hjiv, /Ui)] 72 /Po a 2 (A*i) , 

where L(/x 2 ,//i) = a s (fi 2 )/a s (^i), Po = 11 - |iV>, and 7^ 0) = 50/9. 

The analogous decomposition for the quark-antiquark-gluon matrix element in terms 
of the pion three-particle DA reads : 



(A2) 



(A3) 



X 



(q-x) 



+ 



e^8 P (lxU{(q P 9 5X ~ gV A )^±(«,) + 



In the above 



(ps n = 360aia 2 a 3 



qx(q s x p — q p x 5 ) 
(q-x) 

"l + ^(7a 3 -3) 



(y\\(ai) +(p±(aij) } 



(if\\(ai) +<p±{aa)) } 



(A4) 



(A5) 



is the twist-3 quark-antiquark-gluon DA, with nonperturbative parameters f 3n and de- 
fined via matrix elements of the following local operators: 



(0|wa^7 5 G Q/ 3<i|7r (q)) = if 3v (q 



(A6) 
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3 3 
(0\ua^xl5[Di3,G aX }d - -dpua^xl^Gaxd^' (q)) = —h^^q^q,, ■ 

The scale dependence of the twist 3 parameters is given by: 

y U 7r (/i 2 ) = [L(ji 2 , fil)] ^^{Hl), hnM = [L{^2^l)Y° V 3 ' /37r(/il) 

The corresponding expressions for the twist-3 quark-antiquark DA are: 
<p p (u) = 1 + 30^-Cl /2 (u -u)- 3^f^Cl /2 (u - u), 
Va {u) = 6u(l - «) fl + 5-^- (l - ^) Cl'\u - u)\ . 



(A7) 

(A8) 
(A9) 



(A10) 



_ + 26.(1 - 2a 3 ) 



Finally, the expressions for the four twist-4 DA entering the decomposition (|A4|I are: 
<p\\((Xi) = 120^e 7r («i - a 2 )aia 2 a 3} 
(p±(ati) = 30^(ai - a 2 )a\ 

<p\\{(Xi) = -120<^aia 2 Q!3 ^ + £^(1 - 3a 3 ) 

^(ai) = 305^(1 - « 3 ) + 26.(1 - 2a 3 )l , (AH) 

where the nonperturbative parameters 8%. and e n are defined as 

(0\uG a ^ a d\-K-(q)) = -i8lUq„ , (A12) 
and (up to twist 5 corrections): 



(0\u[D^, G^rfd - -d^uG^^-jx ( g )) = -—fj^q^ ; 
with the scale-dependence: 

SC F 10 



(A13) 



(A14) 
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Figure 3: The ratio of the (7T + 7T~\Os g \B®) and (7r + 7r~|Oi|.B°) matrix elements calculated 
from LCSR (solid line) as a function of the Borel parameters Mi in the pion channel (a) 
and M 2 in the B channel (b). Long-dashed, short-dashed and dash-dotted lines are the 
hard-gluon, soft-gluon and quark-condensate contributions to the sum rule, respectively. 
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